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Abstract 

In this paper, we consider a second order nonlinear ordinary differential equation of the form 
x + fei-^ + (&2 + k^x)x + k^x 3 + k§x 2 + k§x = 0, where k^s, i = 1, 2, 6, are arbitrary parameters. 
By using the modified Prelle-Singer procedure, we identify five new integrable cases in this equation 
besides two known integrable cases, namely (i) £?2 = 0, = and (ii) k\ = 0, &2 = 0, k§ = 0. 
Among these five, four equations admit time dependent first integrals and the remaining one 
admits time independent first integral. From the time independent first integral, nonstandard 
Hamiltonian structure is deduced thereby proving the Liouville sense of integrability. In the case 
of time dependent integrals, we either explicitly integrate the system or transform to a time- 
independent case and deduce the underlying Hamiltonian structure. We also demonstrate that 
the above second order ordinary differential equation is intimately related to the two-dimensional 
Lotka- Volterra (LV) system. From the integrable parameters of above nonlinear equation and all 
the known integrable cases of the latter can be deduced thereby. 
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I. INTRODUCTION AND STATEMENT OF RESULTS 



A. Introduction 

In this paper we investigate the integrability properties associated with the second order 
nonlinear ordinary differential equation (ODE) of the form 

x + ki h (k 2 + k 3 x)x + fc 4 x 3 + k 5 x 2 + k 6 x = 0, (1) 

x 

where fcj's, i = 1,2, ...,6, are arbitrary parameters and the over dot denotes differentiation 
with respect to l t\ Equation (TjQ) is a combination of two different classes of equations, namely 
Lienard type equation {k\ = 0) and equation with quadratic friction {k 2 — k% = 0). Using 
the transformation z = x one can reduce Eq. ([1]) to the Abel equation of the second kind 1,2 . 
However, the resultant equation is not integrable in general. A set of integrable parametric 
choices of this equation has been listed in Refs. 1 & 2. The importance of the study of Eq. 
(CQ) arises from the fact that this equation contains many physically interesting equations 
such as the modified Emden equation, unforced Duffing oscillator, Helmholtz oscillator, etc 
and is intimately related to the much analyzed biological model, namely the two dimensional 
Lotka-Volterra system 3,4 (LV), 

x = x(ai + bux + b 12 y), (2a) 
V = y(a 2 + b 2 \x + b 22 y)- (2b) 

In order to study the integrable properties of Eq. we employ the modified Prelle-Singer 
procedure (PS) which has been intensively used to identify new integrable cases for several 
equations 5,6 . Using this procedure we identify several new integrable parametric choices and 
deduce their corresponding integrals of motion. In order to prove the integrability we either 
associate a conservative Hamiltonian structure to the equation constructed from the time 
independent integrals of motion, thereby proving Liouville sense of integrability, or explicitly 
integrate the integrals of motion to obtain general solution, proving complete integrability. 

During the past three decades or so several attempts have been made to explore the 
integrable cases in the LV system because of the immense importance of the problem in 
mathematical ecology and biology 7-12 . Integrals of motion of the LV system have been con- 
structed for several parametric choices 9 ' 10,12-19 . However, we find that all these parametric 
choices effectively reduce to one of the three general parametric choices for which integrals 
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of motion have been reported 9 . Interestingly, we recover all these three integrable cases in 
addition to many subcases by comparing and expressing the parameters appearing in ([1]) in 
terms of the LV system parameters. The results exactly coincide with the reported ones in 
the literature. Thus, as a by-product, we recover all the known integrable cases of LV by 
investigating the integrability properties of Eq. (pQ). At this point we stress the fact that 
as far as integrability is concerned one can extract complete information about (J2]) from ([1]) 
and not vice-versa. In other words when one tries to deduce the integrable cases of ([1]) from 
02]) only partial results can be extracted. That is one will loose a major portion of the results 
when we proceed in the reverse way (see Sec. VI for more details). 



B. Results 



(4) 
(5) 



Our main results can be summarized as follows. Using the modified PS method we identify 
seven integrable parametric choices of which five seem to be new as far as our knowledge 
goes. The underlying forms of the equations are 

•■ , x 2 n i \ • 7 3 k 2 k±(3 + 2ki) 2 hkK2 + kx) 

x + k 1 — + (k 2 + k 3 x)x + hx 3 + -44q r^ x + J rvk — rr~ x = °> 3 

x k 3 (l + ki) ki(l + kx) 

x + h— + (k 2 + k 3 x)x + ' 3 x 3 + 3 x 2 + k 6 x = 0, 

x (3-\-2ki) z (3 + 2/ci) 

, x 2 , , k 3 k 2 (l + ki) 2 (2 + fci)fc? 

•■ , x 2 , . , 2 2(3 + 2A;i)fc| 

x + k 1 — + k 2 x + k 5 x 2 + ) - 2 x = 0, 6 
x (5 + 4A; 1 ) i 

, x 2 , „ k 2 k 3 9 (2 + k-\)kl , , 

x + h— + (k 2 + k 3 x)x + hx 3 + x 2 + . 2 x = 0. 7 

x [3 + 2kx) (3 + 2ki) z 

For these five cases we prove integrability either in the Liouville sense by constructing the 
time independent Hamiltonian structure or deduce the solution by explicitly integrating the 
time dependent integrals of motion. 

The remaining two cases which are already known in the literature are 

x 2 

x + ki h k^x 3 + k 5 x 2 + k 6 x = (8) 



x 

X 

and 



x + k 3 xx + k A x 3 + k e x = . (9) 
We note that equation (jSJ) can be reduced to the Bernoulli's equation, which can then be 
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integrated to give the solution in terms of quadratures. Equation ([9} is the modified Emden 
type equation with linear forcing term which has been studied in some detail in Refs. 20 
& 21. The Hamiltonian structure for this equation for the parametric choices k 2 > 8A;4 
has been given in Ref. 20. Here we present the Hamiltonian structure for the parametric 
choice k 2 < 8&4 as well. In addition to these seven integrable parameteric choices we find 
the following equation 

x + k 1 — + (k 2 + k 3 x)x+ 3 ; 2 - . J x 2 + k 6 x = o, 10 

x 2(2 + ki) 

where 00 = \Jk\ — 4(1 + k\)ke, for which we obtain a time-dependent integral of motion. 
However, we are able to prove its complete integrability only with an additional parametric 
restriction k 5 = ^0^r^ , and k 6 = t^^jI, which reduces Eq. (TT0|) to Eq. ([5]). In this 
parametric choice we are able to explicitly integrate the time-dependent integral of motion 
to find the general solution. 

The plan of the paper is as follows. In the following section we briefly describe the 
extended Prelle-Singer procedure 5 ' 6 applicable to second-order ODEs. In Sec. IHII we identify 
the integrable parametric choices of Eq. (TjQ) which admit time independent integral of 
motion through the extended PS procedure. In Sec. IV we construct explicit conservative 
Hamiltonians from the time independent integral of motion. Further we transform these 
Hamiltonians to simpler forms using canonical transformations in order to explicitly integrate 
the canonical equations of motion. In Sec. |Vl we then identify the integrable cases of 
([T]) which admit explicit time-dependent integrals of motion. To establish the complete 
integrability of these cases we transform the time-dependent integrals of motion into time- 
independent integrals of motion and integrate the latter and derive the general solution. In 
Sec. |VT] we present the connection between Eq. (CQ) and the LV equation. In Sec. VII we 
rewrite the results in terms of the LV parameters and point out the integrable equations. 
Finally, we present our conclusions in Sec. IVIIII 

II. EXTENDED PRELLE-SINGER (PS) PROCEDURE 

Let us rewrite Eq. ([1]) in the form 

x = —(ki h (k 2 + k^x)x + k^x 3 + k$x 2 + k 6 x) = 4>(x,x). (11) 
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Further, we assume that Eq. (TTTT) admits a first integral I(t, x, x) = C, with C constant on 
the solutions, so that the total differential becomes 

dl = I t dt + I x dx + I x dx = 0, (12) 

where subscript denotes partial differentiation with respect to that variable. Rewriting 
Eq. (fTTj) in the form (pdt — dx = and adding a null term S(t, x, x)xdt — S(t, x, x)dx to the 
latter, we obtain that on the solutions the 1-form 

((p + xS)dt — Sdx — dx — 0. (13) 

Hence, on the solutions, the 1-forms (|12[) and f[T3T) must be proportional. Multiplying ffTBl by 
the function R(t, x, x) which acts as the integrating factor for (1131) . we have on the solutions 
that 

dl = R((f) + Sx)dt - RSdx - Rdx = 0. (14) 

Comparing Eq. (|T2|) with (fl4"]l . we have the relations 

I t = B{<j> + xS), I X = -RS, h = -R. (15) 

Then the compatibility conditions, I tx = Ixt, h± = I±t, Ix± = I±x, between the different 
equations of (fT5j) . provide us the relations 

S t + xS x + (pS ± = -4> x + (p x S + S 2 , (16) 
Rt + xR x + <PR X = -{<p x + S)R, (17) 
R x — SR X — RS X = 0. (18) 

Solving Eqs. (|T6l) - (|T8l) one can obtain expressions for S and R. It may be noted that two 
sets of independent special solutions (S, R) are sufficient for our purpose. Once these forms 
are determined the integral of motion I(t,x, x) can be deduced from the expression 



I = n — r 2 

where 



R + Tx {ri ~ r2) 



dx, (19) 



n 



J R(<f> + xS)dt, r 2 = J(RS+—r x )dx. 
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Equation ([19]) can be derived straightforwardly by integrating Eq. (fl5|) . We solve Eq. 
dl]) through the extended PS procedure in the following way. For the given second-order 
ODE ([T]), the first integral / should be either a time- independent or time dependent one. 
In the former case, it is a conservative system and we have I t — and in the latter case we 
have ^ 0. So, let us first consider the case I t = and determine the null forms and the 
corresponding integrating factors, and from these we construct the integrals of motion and 
then we extend the analysis to the case / t ^ 0. 



III. TIME INDEPENDENT INTEGRALS : INTEGRABLE PARAMETRIC 
CHOICES 

In this section, we identify a set of parametric choices of (EE]) for which time independent 
integrals exist. For this purpose we first find the null forms and integrating factors corre- 
sponding to equation ([T]) and using these functions we construct time independent integrals 
of ([!]) through the relation (fT9l) . 



A. Null forms and integrating factors 

Since I t = 0, one can easily fix the null form S from the first equation in ffl5l) as 

_ -0 _ (fci^ + (fc 2 + hx)x + k A x 3 + k 5 x 2 + k 6 x) 

X X 

Of course one can easily check that the S form given above satisfies Eq. (1161) . 
Substituting this form of S, given in f[2"Ul) . into f[T7|) we get 

±Rx _ ( kl - + (k 2 + k 3 x)x + hx 3 + k 5 x 2 + kex)Ri = ( — - fc ^ 3 + %r 2 + %r \ 

X \ X X J 

Since we are interested in time independent integrals we take R t = 0. As we noted earlier 
any particular solution of the above equation along with the null form S is sufficient to 
construct an integral of motion. To derive a particular solution of (12 lj) we make an ansatz 
for R of the form 

R = {A(x) + B(x)x + C{x)x 2 Y' (22) 
where r is a constant and A(x), B(x) and C(x) are arbitrary functions of their argument. The 
reason for choosing the above form of ansatz is as follows. To deduce the time independent 



first integral I we assume a rational form for 



/, that is, / = f( x ' x \ where / and g are 



arbitrary functions of x and x. Using (|T2|) we obtain S = I x /I± = (f x 9 ~ f9x)/(fx9 ~~ fg±)- 
and R — I± = (f±g — fg±)/g 2 and from these two expressions we find that the numerator of 
R should be the denominator of S and so we fixed the numerator of the R in Eq. ( 1221) as x 
(see expression (1201) for S), that is R = h fe x ) > wnere h * s an arbitrary function. However, it is 
difficult to proceed with this choice of h. So, we further assume that h(x, x) is a polynomial 
in i. To begin with we consider the case in which h is a quadratic function in x, that 
is h = A(x) + B(x)x + C(x)x 2 . Since R is in rational form while taking differentiation or 
integration the form of the denominator remains the same but the power of the denominator 
increases or decreases by a unit order from that of the initial one. So, instead of considering 
h to be of the form h = A(x) + B(x)x + C(x)x 2 , we consider a more general form h = 
(A(x) + B(x)x + C(x)x 2 ) r , where r is a constant to be determined. We note here that this 
form of ansatz played a crucial role in deducing the time independent integrals of certain 
dissipative systems, see for example Ref. 20. 

Substituting ( 1221) into ( 12T1) and solving the resultant equation, we find that the solution 
exists only for certain specific choices of fcj's. In the following we provide these parametric 
restrictions with the resultant forms of R: 



Case (i) : 

k b = ^(3 + 2k 1 )p 1 , k 6 = fc 4 (l + h)(2 + h)pl, pi 



fc 3 (l + fci) 



kx^-l 



x 



((2 + kx)x + A; 3 a;(pi(2 + k x ) + x))x + fc 4 a; 2 (pi(2 + k x ) + x) 



2 ■ 



< 4fc 4 (2 + ki) 



R 



= < 



(2-r)fci 



(23) 



k [k 3 (r - l)(pi(2 + ki) + x)x + (2 + fe^ri] 



k\ > Ak 4 (2 + ki) 
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hi ± hJkl - Ak A {2 + h) 
where r = 



2h(2 + ki 

h 

Case (ii) : 



For the choice ki = — 1, we find k A = 0, k e = — — , R = ; — - (24) 

k 3 k 3 xx + k 5 x 2 



k 2 = 0,k 3 = 0, R = xx 2kl (25) 
Case (iii) : 

fci = 0, k 2 = 0, k 5 = 0, 



x 



. 2k A x 2 + k 3 x(k A x 2 + k e ) + (k A x 2 + k e ) 2 ' 

(x+^k 3 x 2 + ^y 



H < 8k 4 
k 2 > 8fc 4 , r = 



The question now remains as to whether the functions S and R given above satisfy the 
third of the determining equations (vide Eq. ( fTBI) ) or not. One can easily check that all the 
above three sets of functions S and R do indeed satisfy the equation fll8p . The parametric 
restrictions given above fix the equation of motion ([T]) to the following specific forms: 

Case (i) Eq. (TJ): 

x + k-i h fo(pi(l + h) + x)x + k A x 3 + fc 4 pi(3 + 2k 1 )x 2 + k A (l + h)(2 + k x )p\x = 0, 

where p\ 



k 3 (l + h) 
Case (ii) Eq. (TJj: 



x 2 



x + k\ h k A x + k 5 x + k 6 x = 0. 



x 



Case (iii) Eg. (TJj.- 

x + A; 3 xx + k A x 3 + k e x = 0. 

In the next section, we construct time independent integrals for the above equations. 
We mention here that the last equation (Q is nothing but the modified Emden equation 
with additional linear force term which admits a conservative Hamiltonian description for 
all values of k 3 , k A and k e . For the special choice, k A = this equation has been shown 
to exhibit unusual phenomena like amplitude independent frequency of oscillations and 
conservative Hamiltonian structure (for more details one may refer 21). 
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B. Integrals of motion 



Having determined the explicit forms of S and R, one can proceed to construct the inte- 
grals of motion using the expression f[T$j) for the above cases. Substituting the corresponding 
forms of S and R into the general form of the integral of motion ( II 9p and evaluating the 
resultant integrals, we obtain the following time independent integrals for the above three 
cases: 



Case (i): 



log [x 2kl {(2 + h)x 2 + [(2 + k x )px + x][k 3 xx 



+^((2 + fe)» + x)]>]-^g-g 

x kl( - 2 -^[k 3 x{{2 + ki)p! + x) + (2 + h)rx} 
{k 3 (r - l)x((2 + h) Pl + x) + (2 + fc 1 )rx} r_1 ' 

log 



-1 \ ®(x,±) 1 

nfo) J 
PI 



k 3 x[{2 + k 1 )p 1 +x] + 2(2 + kjij x kl 

2(2 + kt)x 
~k 3 x({2 + h) Pl +x) + 2(2 + fe^i' 

^ fc 2 f + fcAfcx + (&| + log[z] - ^ log 
where x) = k 3 x 2 + (2 + k\){k 3 xp\ + 2x), 
and Q(x) = ^4(2 + - fcf (2 + + x x. 
Case (ii): 



Ii 



x 2 



fc 6 2(2 + ki)k b x + (3 + 2^)^ 



x 2 1 + ki 



(2 + fci)(2A;i + 3) 



x 



2(l+fci) 



Case (iii): 



log[2£; 4 x 2 + k 3 (k 4 x 2 + £; 6 )x + (£; 4 x 2 + /c 6 ) 2 ] 



2&3 tan 1 



+ 



x-y^ 8/C4 — fc| 



h = { 



2r 



,f^ 3 X 2 + f ) ^(x+fx 2 + ^ fc3 



( ^2 , r 2 fc 6 



(r-l)fca / [' 



4fc 3 X 



k fcfx 2 + 4fc 3 x + 8fc 6 



logf/cgV + 4fc 3 x + 8/c 6 



A; 2 < 4£; 4 (2 + jfcj) 
k 2 > Ak A (2 + fei) 



(27) 



&f = 4£; 4 (2 + k x ) 
ki = -1 



(28) 



k 3 < 8k 4 



(29) 



fcf > 8k 4 
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Note that in the above the ranges of x and x should be so restricted that no multivaluedness 
occurs. In a recent paper 22 we have shown that the general equation of the form 

x -\ — x + a— — x + X— — - / f(x)dx = 0, (30) 

g(x) g(x) g{xy J 

is related to the damped harmonic oscillator equation 

V" + ay' + Xy = 0, = 4-) ( 31 ) 



dr 



through the nonlocal connection 



y = f f(x)dx, dr = 44^- (32) 
J 9{ x ) 

By suitably choosing f(x) and g(x) one can show that the damped harmonic oscillator 
gets transformed to Eqs. (JHJ), ® and (jSJ). By applying this transformation to the time 
independent integrals of motion of the damped harmonic oscillator one can also obtain the 
above results. However, it is not possible to obtain the time dependent integrals of motion 
of Eq. ([1]) using the above form of nonlocal transformation. We mention here that the 
time independent integrals for the specific equation Q) (Lienard type equation) can also be 
deduced by the procedure described in Ref. 23. 



IV. TIME INDEPENDENT INTEGRALS OF MOTION: INTEGRABILITY & 
GENERAL SOLUTION 

A. Hamiltonian description and integrability 

In the previous subsection we showed that the equations (jHJ) and (jHJ) admit time 
independent integrals of motion. Interestingly one can interpret these integrals as time 
independent (but nonstandard) Hamiltonians for the respective systems and they can be 
treated as conservative systems. In the following we deduce the underlying Hamiltonian 
structure for the Eqs. ©-(El) from the first integrals ( |27j) - (l29l) . To do so we assume a 
Hamiltonian of the form 

H(x,p) — I(x,x) — px — L(x,x), (33) 

where L(x,x) is the Lagrangian and p is the canonically conjugate momentum. From fl33|) 

we get 

91 dp , 

dx dx ' 
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from which we identify 



P 



rdx. 



x 



(35) 



It is clear from ( 1350 that once / is known p can be determined in terms of x and x and 
inverting one can express x in terms of x and p. Substituting the expression for x in terms 
of p into the expression for / and from ( J33l) one can deduce L. Once p and L are known the 
same expression ( 1331) can be utilized to derive H. Using this procedure one can deduce the 
Hamiltonian from the first integrals for all the three Cases (i)-(iii). The Hamiltonians and 
the corresponding canonical conjugate momenta read as follows: 



Case (i): 



log x kl (pi (2 + kx)x + x 2 ) sec 

-^(Piil + kjx + x 2 ), 
where Q(x) = a/4(2 + ki)k 4 - kl 



VLp 



2(2 + fcx) J 



k 2 < 4k 4 (2 + ki) 



(2 + ki) Pl + x 



x. 



(^Jl (par*) £3 _ (pi(2 + kl)x + x 2 ); k 2 > 4A , 4(2 + h) 



(r -2 
log 



P 



r(2 + h) 

+ W?k7) ({2 + h) " l + xh 



kl = Ak A {2 + ki) 



where 



tan" 1 



k 3 x((2+k 1 )p 1 +x)+2{2+k 1 )x 

n 



p = < 



x (2-r)fci 



( i +^Srr(( 2 + fc i)w^ 2 )) 



r — 1 5 



fc| < 4k 4 (2 + ki 



kl > 4fc 4 (2 + ki 



(36) 



Case (ii): 



rr _ r 2( kl+ l) ( P 2 | ^6 x(2(^ + 2)fc 5 + (2fcx + 3)M \ _ ^ 

^ _ X A„2(2k,+D + 7. , 1 H 7-2 , ~I , c },P-*X X.[6I) 



y 4x 2(2k 1 + l) ' ki + 1 



2k[ + 7fci + 6 
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Case (iii): 



H 



( r — 1) (r-2) ( r —\\ , 

(r — 2) 2r 

h ( x 2 4k 6 \ 
| log[p]-yp(^y + -^-J, 




fc§ = 8A; 4 



(38) 



where 



I 



2 



tan 



-i 



k 2 < 8k. 



: 4 , 7 = 7^8£;4 - k 2 



(39) 



X + 



(r-1) 




> 8fc 4 . 



2y 



One can check that in all the above cases the second order equivalence of the Hamilton's 
equation of motion coincides exactly with the associated equations (j3J), © and (jSJ) for the 
appropriate parametric choices. The existence of the time independent Hamiltonian for 
the Eqs. (J3J), (jl]) and (jUJ) assures us that they are Liouville integrable 4 . However, our 
target is to go beyond this statement and construct the solutions for these equations. On 
the other hand, we find that it is difficult to integrate the Hamilton's equations of motion 
corresponding to the above Hamiltonians except for the second case (vide Eq. (137jl ) whose 
solution can be deduced using the standard methods. In order to deduce the solutions we 
introduce suitable canonical transformations so that the new Hamilton's equations of motion 
can be integrated. However, we are able to construct the explicit solution of Eq. (j3J) only 
for the parametric choice p\ = 0, making the resultant equation equivalent to Eq. (J7j) whose 
solution is constructed in the subsection IV CI Hence we do not discuss the solution of case 
(i) separately here. 

B. General solution 

In this sub-section we discuss the method of finding the general solution for the case (ii) 
and case (iii). 



We find that one can rewrite the first integral (128]) straightforwardly into the following 



Case (ii): 
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quadrature 



t-tr 



dx 



T ^.-2fci _ k 4 (3+2k 1 )x 4 +2k ;> (2+k 1 )x3 _ k 6 x 2 
llX 2fe2 + 7fc 1 +6 1+fci 



(40) 



On the other hand the solution of (jHJ) can also be obtained by transforming it to the 
Bernoulli equation through the transformation x = z(x), where x is the new independent 
variable. The solution of this reduced Bernoulli type equation is again given in terms of 
quadratures 2 . 

Case (iii) a): Parametric choice k\ < 8k^ 

In order to obtain the solution for the under damped case (k 2 < 8/C4) we use the following 
canonical transformation, 



x 



— tan 

k A 



h P 

£4 



p 



kJJ 2 
— — cos 



hp 

fc 4 



(41) 



and transform the Hamiltonian H = ?fp (x 2 + j*) — log (x 2 + sec[^(x 2 + f^)] into 

(42) 



H = loe 



The corresponding Hamilton equations of motion are 



fc 6 2 

— sec 


x [h P 


sec 






K4 


V ki 









u 




2\ I - 1 tan 

fc 4 




P 



P 



h 
4 



7 



tan 



7 



U 



(43) 



From the first equation we can express P in terms of U, that is, P 
and substituting this form of P into ( 142]) we obtain 

\4k 6 + htl 2 ) 



^tan" 1 

«6 



H = lot 



4k d 



sec 



= E. 



(44) 



Rewriting the last equation (14"4"I) for U as 



£/ = 2\ /exp 



4 



cos 



7 



[7 



h 



(45) 



and integrating we arrive at the following quadrature 

dU 

t-tn 



2Je W [E+!fU) cos [$U] - % 



(46) 
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Case (iii) b): Parametric choice k\ > 8k^ 

In order to obtain the solution for the over damped case [k 3 > 8/C4), we use the following 
canonical transformation x = p = H^- and obtain a new Hamiltonian of the form 

*-*±-(^ + ¥)- m 

where <j\ = -^7. The canonical equations in the new variable assume the form 

2 r 12 

fr hP{l-r) • hrU 2U^ 

U ~ 2r ' P - — -^~U- (48) 

From the first relation one can express P in terms of U, that is P = ^cl-v) • Substituting 
the latter in the Hamiltonian we get 

H = E = a 1 U^-^U> + -{f- (49) 
2k 3 k 3 (l-r) 

which in turn leads us to the quadrature of the form 

t - t0= (Ml^l))7 £_ (50) 

Case (iii) c): Parametric choice fc| = 8/c 4 

Finally, now we focus our attention on the critically damped case (&§ = 8^4). Using the 
same canonical transformation x = jj, p = H^- we rewrite the underlying Hamiltonian ( )38l) 

as 

= 21og[l/]-^(^ + lfc 6 tf 2 ). (51) 

The corresponding canonical equations turn out to be 

-k 3 P ■ k 3 k 6 2 , , 



Substituting P = — ^ in the Hamiltonian, we get 



H = E = 2 log[U) - ifcgA^ 2 - -^f) 2 (53) 
4 k 3 



Rearranging we get, 



r = . 8k 3 \og[U}-klk 6 U 2 -4Ek 3 (r . () 



^3 - 2 
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which upon integrating reduces to the following form of quadrature 

t-to = V2f - dU (55) 

J y/%h\og[U\ - kJhlP - 4Ek 3 

Summarizing the results, we find that the nonlinear equations (J3J), (jSJ) and (jUJ) admit 
time independent Hamiltonians for all values of the system parameters and can be classified 
as integrable ones in the Liouville sense. While constructing the solution for these three 
equations we find that the nonlinear system (jHJ) can be transformed into Abel equation 
which in turn can be integrated into a quadrature. Solution of Eq. (JHJ) is given in terms of 
quadrature by applying suitable canonical transformation to its Hamiltonian. By using again 
canonical transformations the Hamiltonian corresponding to Eq. can be transformed and 
the canonical equations corresponding to this transformed Hamiltonian can be integrated 
and the solution is given in terms of quadratures for the parametric choice p\ — (see Sec. 



V. TIME DEPENDENT INTEGRALS (I t + 0) 

In this section, we explore the parametric choices for which (JT]) admits time dependent 
integrals. The underlying procedure is same as that of the time independent integral case 
but involves somewhat lengthy calculations. As a first step in this process we derive the 
null forms and integrating factors corresponding to Eq. (CQ) by solving the determining Eqs. 

(HE}- (HE}. 



A. Null forms and integrating factors 

In the previous section we considered the case I t = 0. As a consequence S turns out to be 
However, in the case I t ^ 0, the function S has to be determined from Eq. (flE]) . that is 



S t + xS x h (k 2 + k 3 x)x + k 4 x + k 5 x + k 6 x S± 



+k 3 x + 3hx 2 + 2k 5 x + k 6 -S + (k 2 + k 3 x)^j + S 2 . (56) 

Since it is too difficult to solve Eq. (|56|) for its general solution, we seek a particular solution 
for S, which is sufficient for our purpose. The time independent integral case clearly indicates 
that S should be in a rational form. To begin with one may consider S = ^jfy, where / 
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and g are arbitrary functions. However, it is difficult to solve Eq. ( 156)) with this form of 
S. So one may assume that / and g are simple polynomials in x with coefficients which are 
arbitrary functions of t and x. We seek a simple rational expression for S in the form 

_ a(x, t) + b(x, t)x + c(x, t)x 2 

~~ d(z,t) + e(a;,t)i; + /(a;,t)i; 2 ' ^ ^ 

where a, 6, c, d, e and / are arbitrary functions of x and t which are to be determined. One 
may also consider a cubic polynomial in x both in the numerator and in the denominator. 
However, the resultant analysis did not yield any new result. So we confine our presentation 
here to the form (157)) only. Substituting (157)) into (156)) and equating the coefficients of 
different powers of x to zero, we get 

—a 2 + a t d — adt + adk 2 — d 2 k$ + adk 3 x — 2d 2 k 5 x — bdk e x + aek 6 x — ?>d 2 k^x 2 

—bdk^x 2 + aek^x 2 — bdk^x 3 + aek^x 3 = 0, 
2 2adk\ 

a x d — 2ab — ad x + b t d — bd t + ate — aet + 2aek 2 — d k 3 — 2dek^ H h 2aek 3 x 

x 

—4dek 5 x — 2cdk s x — 6dek^x 2 — 2cdk 5 x 2 — 2cdk 4 x 3 = 0, 



b x d — b 2 — 2ac — bd x + c t d — cd t + a x e — ae x + b t e — be t — cdk 2 + bek 2 — 2dek 3 — e 2 k^ 

d 2 ki bdk\ 3aeki ,, 2 , 

H — I 1 cak 3 x + bek 3 x — 2e k 5 x — cek 6 x 

iXj i\j JLj 

—2>e 2 kiX 2 — cek$x 2 — cek^x 3 = 0, 

2 2dek\ 2bekx 
c x d — cd x — 2bc + b x e — be x + c t e — ce t — e k 3 H 1 = 0, 

x 2 x 

2 e k\ cek\ , . 

c x e-c ~ce x -\ — H = 0. (58) 

or x 

We obtain four parametric choices for which we are able to find nontrivial particular 
solutions of the above set of coupled partial differential equations. We present the parametric 
choices and the null forms S in Table-I. We find that the null forms S obtained are incidently 
independent of time. 

Next, the above identified forms of S along with their corresponding parametric restric- 
tions are substituted in the determining equation ( ITT)) for R. To solve the resultant equation 
for R we make use of the ansatz 

d = F(t)S d , , 

{A{x) + B{x)x + C{x)x 2 Y' 1 1 
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where Sd is the denominator of S. We demand the above form due to the following reason. 
To deduce the first integral / we assume a rational form for J, that is, I = f^' x,x \ where 
f and g are arbitrary functions of t, x and x. From ( lisp , we know that S = I x /I± = 
{fxg — f9x)/{f±g — fg±) and as we find the deduced S forms are independent of time, the 
numerator and denominator of S should share a common factor which is a function of time 
alone, that is, the numerator and denominator of S should be of the form (f x g — fg x ) = 
F(t)fi(x,x) and (f x g — fg x ) = F(t)gi(x,x), respectively. Moreover, from the relation R = 
I± — (f±9 ~ f9x)/9 2 i we find that the numerator of R should be the denominator of S. 

On solving the resultant equation, we obtain the integrating factor R. We present the 
forms of R along with the S forms in Table 1. Once S and R are determined then one has 
to verify the compatibility of this set (S,R) with the extra constraint Eq. ( fl8l) . Having 
verified the compatibility of (S,R), we substitute S and R into Eq. ffl9|) and construct the 
associated integral of motion. In this way we identify four sets of integrable parametric 
choices. We present the explicit form of these integrals of motion in Sec IV Bl 



B. Integrals of motion and general solution 



{l + ki)k% . k 2 k 3 



Case (i) ki, &2, k%, k& : arbitrary, k^ = — , k 



(3 + 2fci) 2 ' '" 5 3 + 2h 
The parametric restriction given above fixes the equation of motion ([T]) as 

, i 2 , k%(l + ki) □ kok^ 9 , m 

x + h— + (k 2 + k 3 x)x + :\ J x 3 + x 2 + k 6 x = 0. © 

Equation Q reduces to the generalized modified Emden equation for the parametric choice 
k\ = whose integrability has been studied in detail in Ref. 5. The integral of motion 
associated with Eq. (0J for k% ^ turns out to be 

fx + Irror-U + P2X 2 \ 

( ia ) A = e^ ( X , h^-l,u^0 (60) 

(ib) Ii = ^ Wi+te?)+*** y ki = _ x (61) 

2(3 + 2k 1 )x 

(lC) 1 ~ * " (3 + 2h)k 2 x + 2(1 + k^ihx 2 + (3 + 2h)x) ' (62) 



where p 2 = ctt, w = y^l ~ 4 ( 1 + fc i) fc 6- 
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TABLE I: Parametric restrictions, null forms S and integrating factors R of x + k\ ^- + (k 2 + &3x)x + /C4X 3 + /C5X 2 + k$x 



Case 


Parametric restriction 


Form of S 


Form of i? 


(i) 


, _ (i+fci)fcg 

M ~~ (3+2fci) 2 ' 

7„ r _ k 2 k 3 
^5 ~~ 3+2fci • 

?v\ , n.2 ; aiuiiiaiy 


k^x x 
(3+2fei) x 


| ~ , 1 fc 2±" | H x i \ 
I x ^2 x 1 3+2fe 1 x / 

, , — /l,2 /J/1 1 J. XL 


(ii) 


^ ~~ 2(2+fei) ' 

fc 4 = 0. 

ki, k 2 , /%3 arbitrary 




e \{k 2 ±u)t x k 1 


(iii) 


, 2(3+2fci)fc| 
ft 6 (5+4fci) 2 ' 

fc 3 = 0, A; 4 = 0. 
fei, ^2 arbitrary 


4(1 r +^ )fc 2 x 2 + (5+4fci)fc5a; 3 +2(l+2fci)A:2xiH-fci(5+4A:i)i 2 ) 


2(3+2fc 1 )fe 2 t / \ 

e 5 + 4fc i x 2fcl ( 2/c 2 x+ (5 + 4fci)x) 


2fc 2 a: 2 + (5+4A;i)a;a; 



Case 


Parametric restriction 


Form of S 


Form of i? 


(iva) 


U u — k 2 k 3 
^5 — 3 +2 fc! ' 

, (2+fei)fc2 
K6 ~~ (3+2fci) 2 • 

^l) ^2, ^4 arbitrary 


, x , i. _ | (3+2fci)fc 4 x 3 (l+fei)fe 2 
M x -r ft-3-t -r (3 + 2fci)x+fc 2 a; 1 3+2fei 


^/c 2 x + (3 + 2fci)i^ (2 + kx)k\x 2 + (3 + 2ki)k 2 x(k 3 x 2 + 2(2 + ki)x) 
+ (3 + 2fci) 2 (A;4X 4 + ±(k 3 x 2 + (2 + fci)i)) 


(ivb) 


7- — k 2 k 3 
^5 — 3+2A;! ' 

, _ (2+fc!)fc2 
ft ' 6 ~~ (3+2fci) 2 ' 

, (r-l)fc 2 
K4 (2+fci)r 2 ' 

^l) ^2i ^3 arbitrary- 


I. i | | (3+2fc 1 )fe 4 a : 3 


(2-r)(2+fc 1 )fc 2 * 
e 3 + 2fc i (fc 2 a;+(3+2fci)i)x fe l( 2 -'') 

^(2+fci )r ( +x) +fc 3 (r- l)a; 2 ^ 



Rewriting (160]) - (|62|) as first order ODEs we find that the resultant equations are of 
Bernoulli type which can be solved using the standard method 2 . The general solution to (j3J) 
in each of these cases turns out to be 



c 2 i 



(ia) x(t) = < c 2 e 2 



xF 



h 



71 



h) 2 " h ~ 2p 2 



±(ci - c 2 ) =Fc 2 1 r 2 ( 



2w 



2w 



2w' Ji 



(ib) = (7 2 e* + fc 3 e' 



71 



71 



„±LOt 



± c 2 



c 2 ^^^ (63) 



7i 



(64) 



(ic)x(f) 



^(3 + 2^) (-c 3 ) 1+fc i exp 



2 



(7i - t) 1+fc i 



2W(1 + fc 1 )fc 3 r 



2 + fci c 3 (t - h 
T+k~i 2 



-7 2 (3 + 2A; 1 )A; 2 (-C3) T w) \ W 



(65) 



where 7 1 is the hypergeometric function, T is the gamma function, c\ = f^^, c 2 = 



c 3 



k 2 
1+Jfci 



_ _ fc2fc s*+e fc 2'Ji ^ an j j 2 are integration constants. 



Case (ii) &i, fc 2 , &3 : arbitrary, k 5 = ^^-ft^ ' &4 = 0, = \/k\ — 4(1 + &i)&6 
The parametric restriction given above fixes the equation of motion (CQ) as 



^(A^icj) 

x + fa h (k 2 + fc 3 x)a; + 

a; 2(2 + fa) 



x 2 + fax = 0. 



The system (TlTJT) possesses a first integral of the form 



x + p 3 x 2 , fci/ -1, 



(66) 



2(1 + ^) 
where p 3 = (^j- 

For the parametric choice fci = — 1, we obtain the same first integral as we have obtained 
in Case (ib) above. On the other hand for the choice fa = —2, Eq. f )66|) gets reduced to the 
Riccati equation which in turn can be integrated by the standard methods and the solution 
can be obtained of the form 2 

2cie Cl * 



x(t) 



(67) 



h + 27 2 cie 2c i* - 2p 3 e c i* ' 
where c% = \{fa =t u), I\ and 7 2 are the integration constants. 

For other choices of fa one is able to integrate the first integral ( 166]) only by imposing addi 
tional parametric restriction. For example, choosing k { 



(2 + fa)k 2 k 3 k 2 {l + fa) 

and k 5 



(3 + 2fa) 2 ' 



(3 + 2fa) 
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one is able to transform the time dependent integral of motion (166]) into the following time 
independent integral of motion 

k 3 (l + ki) ±thx d . . 

where w and z are the new dependent and independent variables, respectively and are given 
by 

(*+*i)*2* 1+fci (3 + 2A; 1 )e -3 +^T 
w = e 3+ 2fcl a-i+fe^ « = -i tl . (69) 

h 

The additional parametric restrictions fix Eq. ( ITUj) to the specific form 

g + jbi 7 + (jb 3 + ib,x)x+ (3 + 2fei) x +(^^^ = 0, © 



Integrating (168)1 we get, 

dw 



z - z 



T fe 3 (l+fcl) SI 1 ' 



1 + ki 3 + 2fci feio^i " 

2 + fci' ' 2 + fci ' Ji 



'r, (70) 



fcq(l -|- fc]) 

where F is the hypergeometric function 24 and k = . We also mention here that 

2 -|- k\ 

Eq. (jSJ) contains several sub-cases which are already known to be integrable, see for example 
Ref. 1. 

Case (iii) fc 2 , &5 : arbitrary, fc 3 = 0, k 4 = 0, k 6 = 2 f 5 !^ ) 2 2 

The parametric restriction given above fix the equation of motion (CQ) as 

x + h— + p 5 (5 + 4fci)i + fc 5 x 2 + 2(3 + 2fci)p!z = 0, © 
x 

k<i 

where ps = -. The first integral for this equation is 

(5 + 4fci) 

h = x 2k - [x 2 + Ap\x 2 + fgg + 4p 5 xi^ e 2 ( 3+2fcl ^ 5t . (71) 

This time independent integral can be transformed into a time independent integral by 
introducing a transformation of the form 

w = t 1+kle2(1+ki)p5t z = _e 

>/2 Ps 
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where w and z are new dependent and independent variables, respectively. Rewriting the 
integral of motion in terms of the new variables, we obtain 



J x = w ' 2 + 2 V aci+fci) 



1+3) fc B (i + h) 2 w 



(3 + 2fci) 

Eq. (1731) can be integrated further and one can obtain the general solution as 



z - z = w. 



h 



3+2k 1 

I\ — kw 1+fc i 



3+2k 1 

1 + Jfei 1 4 + 3fci kw T +^~^ 



3 + 2ki 2 ' 3 + 2ki 



h 



(73) 



(74) 



(3+2fc t ) , , , , 2 

where F is the hypergeometric function 24 and k = 2 2 < 1+fc i> fe^N ■ 



_ (2+fci)fc 2 
6 ~ (3+2fci) a ' 



Case (iv) ki, k 2 , k 3 : arbitrary, k 5 = , k 

The equation of motion in this case turns out to be 

x + ki h (k 2 + k 3 x)x + k^x 3 + 



.r- , , . , k 2 k 3 2 i /.-.]( 2 • /.-, 



(3 + 2fcj 



-x* + 



(3 + 2A;x) 2 



x = 0. 



The first integral reads 



(iva) Ji = log 



x 2kl (2 + An) (x + p 6 x) 2 + A; 3 x 2 (x + p 6 x) + fc 4 ar 



2/c* tan 



-if 



2(2+k 1 )(x+p 6 x)+k 3 x 2 
\ x 2 yj '4(2+k 1 )k 4 -k 2 



yjA{2 + k l )k i -k 2 

x kl( - 2 ~ r ^ (k 3 x 2 + ( y 2 + k 1 )r(x + xp 6 ))e < - 2+kl ^ 2 ~ r ^ t 

(ivd) I\ — — r 

(k 3 (r - l)x 2 + (2 + ki)r(x + xp 6 )) 

(ivc) h = (2 + kjptf + log [k 3 x 2+kl + 2(2 + h)x kl {x + p 6 x 

2(2 + k 1 )(x + p 6 x) 



+ 2(2 + fci)p 6 t, 4h(2 + ki) > k\ 



4fc 4 (2 + ki) < k\ 



Where P6 = (3^ 



A; 3 x 2 + 2(2 + A; 1 )(x + xp 6 )' 

fc§±fc 3A /fc§-4fc 4 (2+fci) 



4£; 4 (2 + jfej) = jfe 



3- 



[)' ~ ~ 2fc 4 (2+fci) 

In the above forms of Ii, we now introduce the following transformation, 

w = xe P6< , z = , 



(75) 



so that in the new variables the integrals of motion read, 

2k 3 tan -1 

(iva) h = log [w 2kl {{2 + k^w' 2 + [k 3 w 2 w' + hw A }}} - - 
M 2 ~r){k 3 w 2 + (2 + kx)rw'} 



k-iW 2 +2{2+k±)w' 
w 2 yjA(2+k 1 )k i -kl 



V4fc 4 (2 + fci)-fci 



(ivb) h 



w 



{k 3 (r - l)w 2 + (2 + ki)rw'} 



nr-l' 



(ivc) Ji = log 



fc 3 w 2 + 2(2 + A; 1 )u; , )w fcl 



2(2 + fe^u;' 



k 3 w 2 + 2{2 + k l )w r 



(76) 
(77) 
(78) 
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In terms of the new variables w and z Eq. ([7]) reduces to the form 

/2 



w" + ki h k 3 ww' + k^w 3 

w 







(79) 



which is equivalent to Eq. fl3]) with fc 2 = 0. We find that it is very difficult to integrate 
the resultant integrals even after removing the time dependent factors. To establish the 
integrability of Eq. ([7]) we transform these time independent integrals into time independent 
Hamiltonian and thereby establish the Liouville integrability of Eq. flTJ. 

By following the procedure given in Sec. IIV Al we obtain the following Hamiltonian for 
the above cases, namely 

^V 4 ( 2 + h)k A - k 2 z p 



(iva) H = log 



w 2+kl sec 



2(2 + jfej 



(ivb) H 



(r-2 
(ivc) if = log 



(r - 1) , I, \ 7— — y &3P(r - 1) 



r(2 + fci) 



4fc 4 (2 + fcx) < fc 



P 



+ 



k 3 pw 2 
2(2 + A^)' 



4fc 4 (2 + fci) = fcg 



(80) 

(81) 
(82) 



where p is the canonical momenta defined by 



(iva) p 



2(2 + An) 



w 



(ivb,c) p 



2 \/4(2 + ki)k 4 - k 

w (?-r)ki 



tan 



k 3 w 2 + 2(2 + kjw' 
w 2 ^/A(2 + k l )k 4 -k\ 



w' + 



r-1 ' 



4fc 4 (2 + fcx) > A; 2 (83) 
4A; 4 (2 + fci) < k\ (84) 



r(2+fci) 

The existence of the time independent Hamiltonian confirms that the system fl7|) is an 
integrable one. However, in the following we briefly point out the method of integrating the 
underlying Hamilton equations of motion associated with the Hamiltonians (1801) - ( 1821) 



C. General solution 



To derive the general solution, we use the following canonical transformations. 
Case (iv a): Parametric choice k\ < 4/c 4 (2 + k\) 

By introducing the canonical transformation, w = p and p = where U and P are new 
canonical variables, we transform the Hamiltonian ( 1H01 (for the choice fc| < 4/c 4 (2 + k{j) to 
the form 



H = E = 4(2 + ki) log 



sec 



f/V 4 ( 2 + fc l) fc 4 - fc; 



4(2 + fci' 



(85) 
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The underlying canonical equations of motion then become 

u'-. 



-4(2 + h) 2 (, d 



P 



dz 



P' = 2faU 



4(2 + k 



U 



— - 2U (^Ja(2 + h)h- kjj tan 



£/ 2 a/4(2 + ki)k A - k 



4(2 + k x ) 



(86a) 



(86b) 



Eq. ( )86|) can be solved in the following way. Expressing P in terms of U, by using the 
Eq. ( I86al) . and substituting it in ( 185]) we obtain 



E = 4(2 + h) log 



m(uu- 



l\2+ki 



sec 



4(2 + h 



- k,U\ 



(87) 



where [i\ = ^xiTW^J anc ^ ^ ls an arbitrary constant. 



(2+fci) 2 ( 2 + fc 

By splitting [/' and U in OH 



exp 



E+k 3 U 2 
4(2+fci) 2 



£7 I //i sec 



C/ 2 v /4(2+fci)fc 4 -fc i 
4(2+fei) 



l ) 
(2+fel) 



(88) 



and integrating the above expression, we get 

1 Aai + ik 3 4ai(5 + 2ki) + ik. 



z - z 



X 



2 + ibi'8ai(2 + A;i)' 8ai(2 + A;i) 

3+fcj -(E+k 3 U 2 )+2ia 1 U 2 
(2 + fcl)22We 2(2+fcI) 



- ( 2-„ r' 



(fc 3 - 4iai)// 1 



i 

2+fel 



(89) 



where ai = ~ an d F is the hypergeometric function 2 . 

Case (iv b): Parametric choice k 2 > Ak 4 {2 + k{) 

The Hamiltonian (1ST]) (for the choice fcg > 4fc 4 (2 + can be rewritten in terms of the 

as 



new canonical variables, w — % and p- '' 



2 • 



H = 2(2 + k x )rr 12 [p^\j-^) rv2 _ 2 ri2 £;3(r - l)t/ 2 , 



(90) 



where we have defined T\2 



(r-1) 
(r-2) 



. The Hamilton equations of motion corresponding to the 



above Hamiltonian are 

_, 2(2 + k x ) 2 rr 2 2 f P^\ n2 



U' 



P 



P' = 2 ri2 k 3 (r- 1)2U + 



2Jbi(2 + Jb 1 )rr£ 2 / p( 2+k ^ 



U 



JJkx 



(91a) 
(91b) 
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Now we integrate Eq. (19 lap by following the same analogy described in the previous subcase. 
First we rewrite Eq. (19 lap for P and obtain 



P 



U k ^ 2 U' 
2(2 + ki) 2 rr\ 2 



(2 + fcl)ri2-l 



(92) 



H = fi 2 U mi U' m2 + u, 3 U 2 = E, 



(93) 



where 



2(2 + h)rr 12 



u „ _ onai. ( r i\ 

(2(2 + k l yrrf 2 ^ {2+kl)rl2, ^ " ^ ^ 
mi = fcj ((2 + fei)m 3 - 7*12) , m 2 = r 12 (2 + k x )m 3 . 



From (I93|) one can express 

[/' 

Integrating ( |94|) we get 

m 2 C/ /> 2 P mi 

Z- Z 



E - fi 3 U 2 



nil + m 2 1 m i + 3m 2 /x 3 P 2 



2m 2 ' vri2 ' 2m 2 P 



(94) 



(95) 



nil + "^2 V E 

where F is the hypergeo metric function 2 and t is an integration constant. 
Case (iv c): Parametric choice k\ = 4/c 4 (2 + fci) 

We use the same canonical transformation w = ^, p = and rewrite the Hamiltonian 



( I82p with the parametric choice fcf = 4(2 + hi) as 



P = loe 



U k 



pki 



log[P 2 ] + 



4(2 + fci) 



P 2 . 



(96) 



The associated canonical equations of motion now become 



P' 
P' 



■p(2 + *i), 



P 4(2 + A*) 

Rewriting (I97aj) for P = and substituting the latter into (196]) we get 



(97a) 
(97b) 



P = log[P fcl (-4P- 



'\2+Jfci 



4(2 + fci 



= E 



(98) 
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which in turn can be brought to the form 

Now integrating the above equation we get 

EC/ 2+fc i 
2 - z o ~ JEEH 

( -k z U 2 \ 2 + fe i 

4 + 4fc i / \ 

where E = —2 2+fc i exp ( ) and T is the gamma function 24 . 

We summarize the results obtained in this section. Solving the determining equations 
( I16p -(ll8 p we find that the system ([I]) admits time dependent integrals for four parametric 
choices and their respective equations are (Bl,fl5]), (jSJ) and (JTj). For (jlj) explicit solution is 
deduced by integrating the corresponding time dependent integral of motion. Using suitable 
variable transformations, solutions of Eq. (jSJ), Eq. ([6]) and Eq. (JTJ) are found in an implicit 
form. 

VI. CONNECTION WITH 2D LOTKA-VOLTERRA SYSTEM (LV) 

The detailed study made on the integrability of the second order ODE ([1]) in the previous 
sections helps one to identify the dynamics of 2D-LV system which we find to be related to 
the second order ODE (0Q) under appropriate choice of fcjS. The 2D-LV system 

x = x(a\ + bnx + buy), 

y = y(a 2 + b 2 ix + b 22 y), © 

models the the population dynamics of two interacting species 3 . Interestingly Lotka-Volterra 
systems arises in other branches of physics also such as the coupling of waves in laser physics 25 
and the evolution of electrons, ions and neutral species in plasma physics. In hydrodynamics 
they model the convective instability in the Benard problem 26 . Similarly, they appear in 
the interaction of gases in a background host medium 27 . In the theory of partial differential 
equation they can be obtained as a discretized form of the Korteweg-de Vries equation 28 . 

Since the 2D-LV system is a planar dynamical system it is also being thoroughly inves- 
tigated from a mathematical point of view. Due to the multi-faceted importance of the 



1 + 



2 + V4(2 + h) 2 



(100) 
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2D-LV system, several in-depth and independent studies have been made to classify the in- 
tegrable cases 7-19 . Integrals of motion of the 2D-LV system (T5]) have been studied for several 
parametric choices all of which reduces to any one of the following 3 parametric choices or 
subcases thereof : 

a, = a tt U \ b i < (102) 

M&12 - M 

u b n b 22 

021 = 7 , (103) 

012 

a x = a 2 . (104) 

for which integrals of motion have been explicitly deduced 9 . Interestingly, we find that all 
these three cases (see LV 1, LV 8, LV 15 below), in addition to several subcases of the above 
parametric choices, and their associated integrals of motion can be deduced from the results 
of Eq. straightforwardly. 

In the following we show that Eq. (T2]) can be transformed to the form (TjQ) and thus the 
integrable cases of ([I]) can be correlated with the integrable cases of (T5]). 

A. Transformation 

To transform the system (T5]) to the form of Eq. (CTJ, first we rewrite Eq. ( l2al) for the variable 

U as 

y = T-(-- b ^ x - ' 012 + ( 105 ) 

&12 / 

Then we substitute the latter into Eq. (l2b]) and obtain the following equation, 

* - ( 1 + - + ( (26h^ - b n - b 21 )x + (2aM - a 2 )) x + (b 21 b n - x* 
V b 12 J x V. o X2 b 12 J V 012/ 

+ fona2 + o 2 iOi - 2ai&n^ J x 2 + \ a x a 2 - J x = 0. (106) 

V " o l2 J V 12 J 

which is of the same form as ([1]). Now comparing fllQ6j) with Eq. ([1]) we find the parameters 
are connected in the following way 

^22 ^22 ^22 

fci = -(l + — ), k 2 = (2a 1 - a 2 ), k 3 = 2b u - 6 U - 6 2 i, 

0l2 012 012 

k 4 = (6 2 i6ii - ir-blx), k 5 = (6ua 2 + 6 2i ai - 2aib 11 -^-), k e = (a a a 2 - r^-a?)- (107) 

012 012 012 
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One may note that both the Eqs. ([I]) and (I106[) contain six parameters and thus one ends 
up with six relations connecting them. Here we emphasize that one can obtain the results 
of the LV equation from the results of Eq. ([T]) straightforwardly, while it is difficult to 
deduce all the results pertaining to Eq. ([T]) from the known results of LV equation. This is 
illustrated with an example. One of the integrable parametric choices of the LV equation (T2J) 
is a2 = ^ffb^zi^r • Upon substituting the above parametric choice in the relations (11071) and 
solving for ^'s we get, k 5 = ^gffl , h = ^jf+t^ and 6 " = ""^S+S^ ~ where 611 
is a real and arbitrary constant. We find that this b\\ is equivalent to r = ~~2^^y^^ 
of the parametric choice k% > 4(2 + hi). However, we have obtained integrals of motion for 
all the three parametric choices k\ > 4(2 + ki) and k\ < 4(2 + ki) which cannot be obtained 
going back from the results of LV. 



VII. RESULTS IN TERMS OF LOTKA-VOLTERRA EQUATION PARAMETERS 

Now we rewrite the results obtained in Sees. HHI -IVlm terms of the LV parameters using 
the relation (11071) . In total we obtain 16 parametric choices in terms of the LV parameters 
(designated as LV 1 - LV 16) and the results are tabulated in Tables HT1 and IHIl corresponding 
to time independent and time dependent integrals, respectively. It may be noted that out 
of the 16 parametric choices, 15 of them reduce to any one of the three parametric choices 
LV 1, LV 8 and LV 15 (given in the last columns of Tables HT1 and HTT1) . while the remaining 
one (LV 6) is the uncoupled case. Also we note that no new parametric choice is obtained 
other than the ones already reported in the literature 9 . 12 > 13 . 15 > 16 as far as the Lotka-Volterra 
system ([2]) is concerned. 

To deduce the integrable choices LV 1 - LV 16 from the results of Eq. ([I]) one can 
derive the corresponding integrals of motion for each one of the cases from the results of the 
second order equation and the relation (I107p . As we pointed out earlier, the integrability 
of the parametric choices obtained in this analysis have already been established in Ref. 
[9,12,13,15,16]. In the following we illustrate the procedure to deduce the time independent 
integral from the results of the second order equation (TjQ) for the parametric choices LV 1 
and LV 2. Integrals of motion for the remaining cases can be deduced similarly and the 
procedure is straightforward. Therefore we do not present the details here. 
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TABLE II: Parametric cases of Eq. (pQ) possessing time independent integrals of motion in terms of the LV parameters 
Case Parametric restrictions Parametric restrictions in terms of LV parameters Integrable choices of the LV equation 
inEq. © 



(:\ l, _ k 2 k 4 (3+2k 1 ) 
W ft 5 - fc 3 (i +fcl ) 

, _ fc 4 fc|(2+fci) 

Ke ~ fcf(i+fcO • 



(i) (a 2 b 11 (b 12 - b 22 ) + 01(621 - 611)622) = 0, 

(ii) ai = a 2 , 6 2 i = 



jjy- aife22(bll— fe2l) 

LV 2 ai = 02, 621 = 



(ii) A; 2 = 0, k 3 = 



fl2 _2a^ 



0, 611 + 621-^2 =0 



LV 3 a 2 



^,621 = 611 -1 



(iii) fci = 0, k 2 = 0, 

h = 



fe + l = 0, 2aife-a 2 = 0, 
6na 2 + 6 2 iai - 2ai6n|g = 



LV 4 612 = -622, ai = -f , 621 = 
LV 5 612 = -622, a\ = a 2 = 



TABLE III: Parametric cases of Eq. (pQ) possessing time dependent integrals of motion in terms of the LV parameters 



Case Parametric restrictions 
in Eq. © 



Parametric restrictions in terms of LV parameters 



Integrable choices of the LV equation 



(i) 



(i+fci)fc| 

(3+2fci)^ ' 

fc2&3 

3+2fci 



fe2l(fell(fol2-2b22)+fe2lb22) 
&12-2&22 

621(01-02) _ q 
6i2 — 2622 



LV 6 6 2i = 0, 612 ^ 

LV 7 fll = a 2 , bn = -^fc 



(ii) 



fc 5 



0, 

fc 3 (fc2±ai) 
2(2+fci) 



7 , 6^622 _ n 



LV 8 621 



Q2"12 (&11 — &2l)+2ai (612621 -611622 )±Oj (fell fcl2+fol2^21 -2611622) _ q Jjy g 

612—622 



_ ^11^22 
6l2 

a 2 , 611 = 0, LV 10 ai = 0, b n = 



(iii) 



fc 3 
k 6 



0, fe 



, »4 



0. 



2(3+2fci)fc| 
(5+4fci) 2 



_ 621 _ bll = o, bub21 _ 



612 

(ai -2<i2 ) (ai 622+202622-02612 ) 
612—4622 



0, 


LV 


11 


ai 


= 2a 2 , 6n = 0, 621 = 




LV 


12 


«i 


02(612-2622) ^ n u n 
- fe2 , b n - U, 021 - U 




LV 


13 


ai 


= 2a 2 , 611 = 621, 612 = 022 




LV 


14 


ai 


= -CL2, hi = 621, 012 = »22 



(iv) fc 5 

fee 



3+2fci ' 
fc 



k _ (2+fci)fc2 



0. 



621(0-1-02) 
612-2622 

(Ol -02) (ai622+02 (622-612)) _ r, 
6i2-2622 ~~ U 



LV 15 ai = 02 

LV 16 01 = a2(fcl b 2 2 ; b22) , 621 = 



Considering the integrals of motion of the second order system (JT]) we note that case (i) 
have three types of integrals (vide Eqs. fl27j) ) depending on the values of the parameters. 
To rewrite these integrals for the first order LV system §2§ first we check whether the LV 
parametric relations LV 1 and LV 2 are consistent with these conditions. While verifying 
this we find that both the LV systems are subcases of the overdamped parametric choices 
kl > Ak^ + h), see Eq. (j2Tj). 

Once the respective integral has been identified then one can replace the variable x in 
terms x and y (vide eq (T5])) in (1271) . The integrals of motion for the above two LV cases turn 
out to be, respectively, 

^22 (^21 — ^11 ) ^22^11 — ^12 ^21 

LV 1 h = yxh^=^ (ai(b n - b 2l ) + b n (b n - b 21 )x + b u (b 12 - b 22 )y) TO^iil (108) 

f>22 — ^22 

LV 2 I x = (ai + b 22 y){b n x + (b 12 - b 22 )y)^^x^^ (109) 

We also note that the well known LV equation 

x = x(ai — b 2 iy), (110a) 
y = y(a 2 + b 21 x), (110b) 

is a subcase of LV 1. The integral of motion corresponding to (IllOp is deduced from the 
integral ff27|) (with ki = —1) as 

h = b 21 (x + y) + a 2 log[x] - ai log[y]). (Ill) 

For the remaining cases LV 3- LV 16, similar analysis can be performed straightforwardly. 

VIII. CONCLUSION 

In this paper, we have investigated the integrability properties of Eq. ([1]) and shown 
that it admits a set of integrable parametric choices. To identify them we have divided our 
analysis into two categories, that is systems which admit time independent first integrals and 
equations which possess time dependent integrals. After carrying out the detailed analysis 
we found that there exists a new equation which admits time independent integral. To 
interpret this integral as a Hamiltonian we first deduce the corresponding Lagrangian and 
then construct the Hamiltonian using the Legendre transformation. Since we have identified 
a conservative Hamiltonian description for a dissipative system, we expect the study can be 
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extended to the quantum case as well in future. The other two systems which admit time 
independent integrals are already known in the literature. However we have also given the 
Hamiltonian description for both of them. 

We then moved on to identify the systems which admit time dependent integrals. Our 
results show that there exist four integrable cases in ([T]) that admit time dependent integrals. 
We have also reported the explicit forms of these integrals. For the first three equations 
we have also found the general solution from these integrals. Since the integral of the 
fourth equation turned out to be a very complicated one it became difficult to integrate 
it straightforwardly. So first we have transformed the time dependent integral into a time 
independent one. Then from the latter we identified a Hamiltonian. We then introduced 
a canonical transformation to this Hamiltonian and transformed the latter into a relatively 
simpler Hamiltonian. This Hamiltonian has been integrated to obtain the general solution. 

We have transformed the identified integrable choices of the second order equation to the 
LV system. Out of the 16 LV parametric choices obtained, 15 reduces to any one of the 
three parametric choice LV 1, LV 8, LV 15. The 16 th one, namely LV 6, is an uncoupled 
case. Interestingly our results reproduce all the known integrable cases of the LV system in 
the literature. 
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